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Te INTRODUCTION 


Although many thin-wire computer programs have been 
developed for the purpose of analyzing antennas and 
scatterers, few of these programs have been directed toward 
the student of electro-magnetic theory. The majority of the 
programs are directed to the engineer or advanced student 
for the purpose of analyzing designed structures or 


verifying experimental data. 


The purpose of the study is to develop a computer 
program by modifying an existing computer code which can be 
utilized as an educational method to develop insight into 
radiating structures by the beginning student of 


electro-magnetic theory. 


The modified Ohio State University Antennas-Scatterers 
Analysis Program (OSUMOD or ASAP) is directed toward the 
beginning student who does not yet have the expertise | 
necessary to manipulate the input data for proper execution 
of the larger more comprehenSive analysis program. Even 
though ASAP is small in core requirements and is fast in run 
time, it is capable of analyzing structures to assist the 


engineer with design problems. 


Since the resulting program, ASAP, is primarily directed 
toward students, the program has been limited to structures 
which contain less than 50 monopoles (segments), no longer 
than one-fourth of a wavelength, and which have less than 51 
nodes (intersections and endpoints). If a ground plane, 
either perfect or finite is present; the stated limits above 


are halved due to the generation of an image structure. 


IIe ORIGINAL PROGRAM 


A. THEORY 


Reference 1 presents the electro-magnetic theory for the 
analysis of antennas and scatterers in an isotropic, linear, 
and homogeneous ambient medium. The analysis is performed in 
the frequency domain with an excitation caused by either a 


generator or an incident wave. 


In the analysis, a piecewise-sinusoidal expansion is 
used for the current distribution. The matrix equation 
Z I= V is generated by enforcing reaction tests with a set 
of sinusoidal dipoles located in the interior region of the 
wire. Since the current distribution has the same form as 
the expansion mode, this formulation is known as the 


"Sinusoidal reaction technique". 


B. COMPUTER PROGRAM 


Reference 2 presents the computer program corresponding 


to the theory presented in Ref. 1. 


1. Input Format 


In the program, the input data must specify the 
frequency, wire radius, wire conductivity, the parameters of 
the exterior medium, coordinates of the points to describe 
the shape and size of the wire configuration, a list of the 
wire segments, and the indicators for the various outputs. 
Table 1 is the input data necessary to analyze a half-wave 


dipole. 


2. Output Format 


In the original form, the only outputs which could 


be requested by the input data stream are the following: 


ae Antenna Problems 


(1) Current Distribution on the Structure. 
(2) Input Impedance. 

(3) Radiation Efficiency. 

(4) Near-Zone Field. 


(5) Far~Zone Field. 
b. Backscattering Problems 


(1) Absorption Cross Section. 
(2) Scattering Cross Section. 
(3) Extinction Cross Section. 
(4) Complex Elements of the Polarization 
Scattering Matrix 
c. Bistatic Scattering Problems 


Echo Area. 


Table 2 is an example of the output data available for 
data of table 1. 


3. LIMITATION 


Although the program can analyze a structure with up to 
50 segments, 55 points and 60 dipoles modes; it can not 
analyze a structure in the presence of a finite ground 


plane. 
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AN EXAMPLE OF THE INPUT DATA FOR THE ORIGINAL PROGRAM 


TABLE 1 
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AN EXAMPLE OF THE OUTPUT DATA FOR THE ORIGINAL PROGRAM 


TABLE 2 


IIITe MODIFIED COMPUTER PROGRAM 
A. Input Format 


As illustrated in table 1 the format for the input data 
cards is not self explanatory. This format can be determined 
by referring to the FORMAT statements of the program of Ref. 
2. Since the modified program is directed toward the 
student, the input data format was changed to allow free 
format. Reference 2 was written in a form which permitted 
modifications to allow flexibility in specifying input data 
for the analysis program. Appendix B, titled "User's 
Manual", discusses the input data cards necessary for proper 
execution of an analysis problem. Appendix B is self- 
contained and may be used independently of the remainder of | 


this document. 
B. Output Format 


In the original computer program, the absence of labels 
encumbered the output data and lessened the usefulness of 
the program. To improve the usefulness of the modified 
version, detailed labels were added to the output data. As 
With the input data, Ref. 2 was written in a form which 
enabled modification to allow more specific output data for 
the analyzed problem. With the addition of the polar 
plotting package, the far-zone electric field intensity 
polar radiation and reradiation patterns can be plotted. A 
Sample problem can be found on page 120 in Appendix B, 


User's Manual. 
C. Finite Ground 
To enable the student or the engineer to have an 


improved analysis program, the finite ground effects were 


added to ASAP. The theory corresponding to the ground 
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effects, which utilize Fresnel reflection coefficients, is 
discussed in Appendix A, titled "System Manual". Also 
discussed in Appendix A is the modified computer program and 
the corresponding theory. The electro-magnetic theory was 
developed in Refs. 1, 2, and 3; and it is restated with its 
corresponding computer code to assist in the understanding 
of the methods applied. Appendix A is self-contained and 


may be used independently of the remainder of this document. 


IVe CONCLUSION 


The addition of ground effect techniques to the original 
program did not alter the accuracy or the computational 
capabilities of the program. The ground effect techniques 
utilized the results of the original program and modified 
these results to account for the effects of the presence of 
the finite ground. 


To verify the numerical results of ASAP, the input 
impedances of both a horizontal and a vertical dipole were 
compared to the solutions of the exact form of the 
Sommerfield's equation. As can be seen in table 3 the 
finite ground treatment of ASAP agrees favorably with 
Sommerfield's solutions. The ASAP finite ground results are 
also in excellent agreement with the previous computer 


solutions of Refs. 4 and 5. 
Ve RECOMMENDATIONS 
Although the program is a general analysis tool for 
students, several future modifications will enhance the 


program as a design tool for engineers. These items include: 


varying the wire radius on the structure; incorporation of 
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varying 


hon-cadiating elements such as transmission lines 


conductivity, and dielectric 


lation radius, 


insu 


the wire 
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eed and reduce the core requirement 


No attempt was made to utilize the symmetry 


the 


in 
admittance matrix when the ground plane is present. 


Ie 


Symmetry were applied, the structure size limit with the 


ground plane present would be approximately that of the 


without the ground plane. 
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SYSTEM MANUAL 


INTRODUCTION: The Antennas-scatterers Analysis Program 
(ASAP) for thin wire structures in a homogenous conducting 
medium performs a frequency domain analysis of antennas and 
scatterers. The program is applicable in the presence of a 
ground either perfect or finite. This appendix will describe 
the computer program which accomplishes this. Although the 
program waS written for the IBM 360 computer system it can 


be executed on another system with minor modifications. 


A piecewise-sinusoidal expansion is used for the current 
distribution. The matrix equation ZI = V is generated by 
enforcing reaction tests with a set of sinusoidal dipoles 
located in the interior region of the wire. Since the test 
dipoles have the same current distribution as the expansion 
modes, this may be regarded as an application of Galerkin's 
method. Rumsey's reaction concept was most helpful in this 
development, and therefore the formualtion is known as the 


"sinusoidal reaction technigue". 


The main routine and each subroutine is discussed 
separately in this appendix. The writeups for the 
subroutines are arranged alphabetically by subroutine name 
after the main program. Each of the discussions includes the 
purpose of the subroutine, brief description, and a listing. 
After the subroutine writeups is a table of the more common 


symbols used in this progran. 


The input data and program limits are discussed in 
detail in the next appendix titled "USERS MANUAL". 


17 


GROUND EFFECTS: In the modified antenna analysis 
computer program finite and infinite ground effects were 
added by using the reflection coefficient technique. The 
method in which this technique was used required the 
generation of an image structure. In this section the 


reflection tecnniuue will be discussed in detail. 


In order to apply ground effects to the electric field, 
the field for the image structure was first calculated as if 
a ground were not present. Then, the field was decomposed 
into parallel and perpendicular components. (A parallel 
component 1s the component which 1s parallel to the plane of 
incidence. A perpendicular component is one which is 
perpendicular to this plane. The plane of incidence is the 


plane containing the normal to the reflecting surface and 


the incident ray.) 
a 1 
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Consider an image monopole with the electric field in 
the 1 direction. The ray, @, is a vector which is 
perpendicular to 1 and passes thru the point of interest. To 
apply reflection technique, the plane of incident must be 
found. It is advantageous to define a new coordinate system 
(asB,z) where and are parallel to the xy plane with ain 


the plane of incident and B perpendicular. 


If the direction cosines (cosx, cosy, and cosy) are 
known, it can be shown that the components of the field in 


the af (xy) plane have the following relationship: 


E cos $ Sin 6 E 
x 

ue 

E Sin 6 -cos E 
Y 


where @ = arctan (cosy/cosx). 


Now the reflection coefficients for the interface can be 


applied as: 


wtt(R) 
1 
R 
2! a 2 
L 


where R and R will be defined later in this section. 
1] 


Applying the matrix equation above yeilds: 
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R R 
g cos 6 sin 6 oe 


(R) 4 (R) 
E 


E sin 4 -cos 6 
y 


(the square matrix 1S unigue, in that, the inverse is equal 
to the original matrix). Since the image direction is 


opposite to the original nonopole, that is, 


(l xz) .. ee Xe ’ 
original image 


the z component of the field, which is in the plane of 


incident, is given by: 


From electro-magnetic theory the reflection coefficients 
for the fields in medium (1) at the interface with another 


medium (2) are defined as: 


for parallel 


cos 8 - \’ &' - sin2 98 
R er 
H cos 8 + \/f &' - sin2 9 


and for perpendicualr 


i ~ &' cos 6 + \V &!' = sin2 98 


V &*§* cos 98 + &' =~ sin@ Q 


where 8 1s the angle of incident as measured from the normal 


to the interface and 
| = Oo 5 s 
Bt = (6 + F/jO)/(E, + 0 /jo) 


where the subscripts correspond to the mediums above. 
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To determine the relationship between ene R and Be R 


a perfect ground (& = 0, GY=om) was investigated. 
1G 
limit R = -1 
H 


+7 


limit R 
Vv 


But, for a perfect ground the contributions to the field 
from the image monopole would be equal to the field of the 
original monopole but opposite in sign due to the chosen 


reference direction, 


{i (R) {| 1 | 
E = R E = - §£ 
1] 
t (8) iL tL 
E = R E =- § 
al Ie 
therefore 
R = R 
{| H 
R =- R is 
AS Vv 


In summary, the contribution to the electric field of a 


monopole over a ground plane at a given point is given by: 


(R) (R) (R) R 
E = £E cosx + E cosy + £E COSZ 
x y Z 
where 
(R) 
= R E cosx + (R -~ R ) E cosx cos? 6 

X ai i t 

+ (R - R) E cosy sin $ cos @ 


21 


= R E cosy — (k ~ R ) E cosy cos? $6 
V4 i | 11 L 


+ (R - R ) E cosx sin $ cos $6 


_ ny E cosz 


where E is the field without the ground plane present and 


- - &' cos 8 + YY &! = sin@ 8 
1) §©360hrd&*§ cos 0 + YW &* = Sin2 6 


cos €@ - VY &' - sin? 8 
i cos @ + /Y &* = sin 6 


&'§' = & - j(0/& w) 
te 0 


De 


MAIN 


PURPOSE: to control the input, output, and the flow of 
calculations. 
METHOD: The main program controls the flow of the 


required calculations by calling only a few subroutines. 
These subroutines in turn call other subroutines which 
actually do the required calculations. The order of the 
calling sequence is diagramed after the listing for the main 


prograa. 


The DIMENSION statements at the beginning of the main 
routine provides the required storage for a wire structure 
with up to 50 segments, 60 nodes and 60 dipoles without the 
presence of a ground plane. If a ground plane is present 
one-half of the reserved storage is required for the image, 
therefore a wire structure with up to 25 segments and 30 


nodes can be analyzed. 


NM denotes the actual number of monopoles (Segments), 
INM is the corresponding dimension, and the dimension for 
CG,. VG, and ZLD is twice INM. The second subscript for MD 
always has a dimension of 4 to correspond to the number of 


segments meeting at a given node. 


N denotes the number of simultaneous linear equations 
and ICJ is the corresponding dimension. The dimension for C 
is (ICJ * ICJ + ICJ)/2. 


In the statements above statement 4, the initial 
conditions and defaults are established. After calling 
subroutine READ to determine the input parameters, the IF 
statements output the parameters to be used for the 
calculations. In the DO LOOP ending at statement 7, the 


the input data of the structure geometry is stored in order 
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to recall if the structure is to be moved for ground plane 


calculations. 


After the image structure is generated and structure 
location is moved, subroutine SORT is called to determine 
the dipole modes. Prior to calling SGANT, the load and 


generator information is established. 


Subroutine SGANT is then called to calculate the 
elements of the impedance matrix. If FEEDS or GENERATORS 
are specified by the input data stream, subroutine GANT1 is 
called to solve for the current distribution due to these 


forcing functions. 


In the DO LOOP ending with statement 29, subroutine 
GNFLD is called to calculate the near-zone field for the 


current distribution of the subroutine GANT1. 


The subroutine GFFLD is called for the far-zone field of 


the current distribution of the subroutine GANT1 in the DO 
LOOP ending at statement 35. The subroutine GFFLD is called 
again in DO LOOPS ending at statements 42 and 51, if 
bistatic and backscattering calculations are requested by 
the input data strean. 
CAELS (1.0: GANT1 

GFF LD 

GNF LD 

POLPRT 

READ 


SGANT 


SORT 
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BLNK 


PURPOSE: to compress data to the left by removal of 
the blank spaces on the input data cards. 


METHOD: A(I) character is compared to the Dlank; and 
if it is true, the A(It1) character is shifted to the A(I) 


position. 


CALLED BY; READ 
SALLS TOs NONE 

SUNROUTINEG O8LNK ILA) i 

OIMENSLTON A880) 2 

OATA BLANK/® ‘%f/ 3 

K = O 4 

C 5 

DO ft 1=1,80 6 

Joe [-K ? 

eal es A({I) 8 

C 1 tf (ACE). EQ.BLANK) Keke] to 

1€ «{K.EQ.0)} RETURN ll 

A(81-K) = BLANK 12 

RETURN 13 

END 14 
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CBES 


PURPOSE: to calculate the quantity B01 where 


BO1 = J (z) / 3 (z).- 
QO 1 


METHOD: If the absolute value of the argument for the 
Bessel function is less than 12, BO1 1S calculated via the 
power series expansion for the Bessel function in the DO 
LOOP ending at statement 3. If greater than 12, the 
asymptotic expression is utilized at statement 4. If the 
Magnitude of the complex part of the argument for the Bessel 
function is greater than 20, B01 1is set to (0.,-1). If the 
complex part of the argument is negative, the sign of R01 is 


changed prior to returning to the calling progran. 


CALLED BY: SGANT 
CALLS TO: NONE 

SUBROUTINE CBES 12,801) 1 
COMPLEX ARG. CC sCSrEX 2 
COMPLEX BO1,2, FERMJeTERMNgMZ24— UNG2) 3 

OATA PE/3.141859/ 4 

IF (CABS(Z).GE.12-0) GO TO 4 5 
FACTOR = 0.0 é 

TEKMN = 10.,0.) 7 

M224 © -0.2592¢2 B 

TERMS = (1.0,0.0) 9 

€ 10 
N= NP-) 12 
JNONP) © TERMS 3 

he 0 4 

1 M = Mo) 15 
TERMS = TERMI®MZ24/FLOATIM*® NOM) } 16 
JNUNP}) = JNUNPP*TERMD 7 

JF (NP.NE*)) GO TO 2 1B 
FACTOR = FACTO2[¢).O/FLOATIM) 19 

TERM = TERMNe«TERMJ*F ACTOR 20 

2 ERKOR = CABSLTERMJ) 2) 

jf (ERRUR.GT.1-0€-L0) GO TO 1 22 

3 TERM) = 0.582 23 

C 24 
BOL © JNUDI/ INEZ) 25 
REFURN 26 

& Y «© ALMAG(2Z) eT 

YF JABSTY).GT.20.) GO TO $ 2B 

ARG © 1.0,1.)82 29 

Ex = CExPlAaG) 30 

CO = EXO) EX 3) 

cs teal ).0 —~).P@VEX-1L./EX) 2 

BOL = (C$eCCI/ICS-CC) 33 
RETURN 34 

5 Bul = ).0 -).) 35 

LF UYet TsO.) BO) © 1.0eke? 36 
RETURN 37 

ENO 38 
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PSHELL 


PURPOSE: to calculate the mutual impedance tern 
contributed by the dielectric insulation on the surface of a 


thin wire. 


METHOD: The contribution to the impedance matrix is 


calculated utilizing the equation below 


& - &) Iin(b/fa 
Z oe eee ee als 2 a rc LO 
mn 2n jw&s m n 
2 m,n 
where Z is defined in subroutine SGANT, oo is the 


mn 
dielectric constant of the insulation, bis the outer radius 
of the insulation, a is the inner radius, & is dielectric 
constant of the external medium, and F is the Sinusoidal 


expansion function. 


CALLED “BY : SGANT 

CAs 5-1 Oe NONS 
SUBROUTINE OSHELL (AM, B8M,0K,CGOS, SGOS,EPZ2,EP ETA, GAMePI1,P12) e 
COMPLEX CGOS SGDS oe EP2 EP sETAsGAMy PLL P12 900, CST } 
DATA PI/3.14159/ 3 
GO = GAM*OK & 
CST © CEPZ-EP) FE TAPALOG(BM/AM) / (4. 8P[ SEP 2%S$6059S6D05S) 5 
Pll = -CSU*{(GO*#S$GOS*CGOS)} 6 
Pl2 = CSF#(GO*CG0S$+*S$G05$) T 
RETURN 8 
ENO 9 
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EQUAL 


PURPOSE: to determine position (location) of the equal 


Symbol on input data card. 


METHOD: The character search begins in the column 
passed to the subroutine. On returning to the calling 
program, the argument passed is the column following the 


equal symbol. 


CALLED BY: READ 
’ 
CALLS TO: NONE 
SUBROUTINE EQUAL UN) 
COMMON /A/ A(B0) 
DATA EQUAL/#=¢/ 


WR OO Ob OUP 


=e = — 
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EXPJ 


PURPOSE: to calculate the exponential integral with 
complex limits. 





METHOD: The exponential integral is defined as: 
V2 
“Vv 
W12 = & dv = EB (V4) = E (v2) + 420 
V 
Vi 


where the integration path is the straight line from V1 to 
V2 on the complex v plane and 


= 
E (Zz) = | e Gitvenre 
- “ie 
z 


The integration path is a horizontal line in the w plane 
Or an inclined straight line from V1 to V2 the v plane. The 
integer n is zero unless this path intersects the negative 
real v axis at a point between V1 and V2. When there is 


such an intersection, 


je) 
— 
= 

il 


1 if Im(V1) > Im(v2) 


o 
| 
i 


- 1 aif Im(V1) < Im(V2). 


The term j2nn is caluclated below statement 12. 


CALLED BY: GGMM 


CALLS TO: NONE 
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GANT1 


PURPOSE: to consider the wire structure as a 
transmitting antenna and calculate the input impedance and 


current distribution. 


METHOD: Tf a wire antenna is driven by a voltage 
generator v_ located at one of the current sampling points 
i 
1 and if displacement currents are neglected, Ampere's law 
i 
yields 


y =v or, } 
m 1 miu 


where F is the sinusoidal expansion function. Thus, the 


excitation voltages V will vanish everywhere except where 
m 


v is not zero. 
1 


The DO LOOP ending with statement 50 uses the delta-gap 
model defined above to determine the excitation voltage 
CJ({I) for all the dipole modes. These are stored 
temporarily in CG(I). Then subroutine SQROT is called to 
obtain a solution of the simultaneous linear equations. 


SQROT stores the solution (the loop currents) in CJ(I). 


In the DO LOOP ending at statement 80, the complex power 
input and input impedance(s) are calculated. The 
time-average power input (PIN) is the real part of the 


complex power input. 


Subroutine RITE is called to make the transformation 
from the loop currents to the branch currents. If IWR is a 
positive integer, RITE will write out the list of branch 


currents. 


38 


the time-average power 


7 


te GDISS to obtain 
dissipated in the lumped loads and the imperfectly 


MAIN 
GDISS 
RITE 
SQROT 


Finally, GANT1 calcudates the radiation efficiency by 


calling subrouin 
conducting wire. 


CALLED BY 
CALLS TO 
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GDISS 


PURPOSE: to calculate the time-average power 
dissipated in the imperfectly conducting wire and in the 


lumped loads. 


METHOD: The time-average power dissipated by the wire 
is calculated in the DO LOOP ending at statement 1 utilizing 


the equation below: 


R / x 
P = Ga 
d 2na 4 


where R is the surface resistance of the wire and a is the 
S 


radius of the wire. 


The power dissipated by the lumped loads is calculated 
by the DO LOOP ending at statement 3. If the wire is 
perfectly conducting, CMM <0, the first calculation is 


by-passed. 
CALLED BY: GANT1 
CALLS TO: NONE 
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GFF 


PURPOSE: to calculate the far-zone field of a 


Sinusoidal electric monopole. 


METHOD: If an electric line source has length d and 
endpoints at (came 47 2) and Vac eer then the coordinates 


of any point on the source are 


x = x + 1 cosx 
y = y. + 1 cosy 
Z= y +1 cosz 


where cosSx, cosy, cosz are the direction cosines of the l 
axis, and 1 is the distance along the source measured from 


the endpoint (XY ye2)- Let the current distribution on 
the monopole be 


I sinh y(a - 1) + I, sinh yl 


I(1l) = ante 
@) Sinh yd 





where ~ and I are the endpoint currents. The far-Zone 


field of this source is 





a = (cosx cos®@ cos$ - cosy cos@ sing + cosz sin@) a 
- = (-cosx sin@ + cosy cos®@) a 
where 
ne Ygd Ytve2 
E = [ (e - g Sinh Yd - cosh yd) Ie 


1 4nr (1-g2) sinh yd 


ya | Weta) 
+ (e + g sinh Yd - cosh j)yd) I, 2 J 
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£¢1) = x Sin@ cos + im sin@ sing + z cosé 


£2) = x sin® coso + 5 sin@® sind + AS cosé 
g = cosx sin8 coso + cosy sin@® sing + cosz cos®@ 


and (r,8,%) are the spherical coordinates of the observation 


point. 


In this subroutine the range dependence has been 
suppressed. The far field vanishes in the endfire direction 


where GK = 0. If a ground plane 1s present (IGRD > 0) the a 


equation above is decomposed into the x, y, and z components 


and the reflection coefficients are applied before oe and E 


field components are returned to the calling progran. 


CALLED BY: GFFLD 


CALLS EG. NONE 
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GFFLD 


PURPOSE: to calculate the far-field for the thin wire 
structure. 
METHOD: The far-field for the structure is calculated 


from the loop currents. The loop currents are either the 
currents produced by the transmiiting antenna calculations 
of subroutine GANT1 or the currents produced by an incident 


plane wave. 

If the incident field is generated by a distance source 
with spherical coordinates BSA OU a) G the excitation 
voltages induced by a incident plane wave are 


V = [2 dl 
fe ml 


et 


where 


| 
tl 


E ex ae dl 
0 P(Y 0 ) 


where = is a vector constant, ar 1s a vector from the 


coordinate origin to the distance source, and fF is the 


radial vector from the origin to the observation point. 


The field - is generated by test dipole m when 
radiating in the homogeneous medium. USing the vector 


potential, the field at the distance point cae eae is 


-Yr 


jwue 0 ae 
E = —aemact exp(Y r fF) dl 
m 


nm Qnr 
0 


m 
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where the radial component is to be suppressed. From the 


above equations, 


Qnr Yr 
V=ze-\7— 2 e See. 
m jwu o an 


If an antenna gain calculation is desired, INC is set to 
zero. PH and TH denote the spherical coordinate direction 
of the distance observation point. The phi-polarized (EPPS) 
and the theta-polarized (ETTS) components of the electric 


field intensity are returned to the calling program. 


If INC = 1, a backscattering calculation is desired. In 
this case PH and TH denotes the incident angles for the 
incident plane wave. These are also the spherical 
coordinates of the distance source. The outputs returned to 
the calling program include absorption, extinction, and 
scattering cross section for each polarization; scattered 


electric field; and echo areas. 


If INC = 2, a bistatic calculation is desired. In this 
case PH and TH denote the spherical coordinate of a distance 
observer. Since this calculation uses the induced loop 
currents (EP and ET), a backscattering call must preceed 
this calculation. The outputs returned to the calling 
program consist of the scattered electric field components 


and echo areas. 


EPP(I) and ETT (I) denote the phi-polarized and 
theta-polarized far-zone fields of dipole node I with unit 
terminal current. In a backscattering Situation, the 
excitation voltages EP(I) and ET(I) are obtained by 
multiplying EPP and ETT by the constant CJI. Then calls are 
made to SQROT which stores the solution (the induced loop 
currents) in EP(I) and ET(I). RITE is called for the branch 
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currents CG(J), and GDISS is called for the time-average 
power dissipated in the imperfectly conducting wire and the 
lumped loads. This power is denoted PDISS and TDISS for 
phi-polarized and theta-polarized incident waves, 


respectively. 


In scattering problems, the incident plane wave has unit 
electric field intensity at the origin. GGG denotes the 
time-average power density of the incident wave at the 
origin. ACSP and ACST denote the absorption cross sections 


for the phi and theta polarizations. 


PIN and TIN denote the time~average power input to the 
Wire structure, delivered by the equivalent voltage 
generators VP and VT at the terminals. PIN and TIN apply for 
the phi and theta polarizations, respectivity. The 
time-average power input is regarded as the sum of the 
time-average power dissipated and the time-~average power 
radiated or scattered by the wire. ECSP and ECST denote the 
extinction cross sections and SCSP and SCST denote the 


scattering cross sections. 


The distance field is calculated in the DO LOOP ending 
with statement 7 for scattering situations, and in the DO 


LOOP ending with statement 9 for the antenna situation. 


The radar cross sections (echo areas) SPPM, SPTM, STPM, 


and STTM, are defined as 


4: Zee 
o= limit 4nr e Sr (Ss 
r S 1 


where S and S_ denote the time-average power densities in 
Ss 1 


the scattered and incident fields evaluated at the origin. 


For an antenna, the following definition is employed for 
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the power gains: 


a. Zeca 
G (8,¢) = limit 4nr e Sr, O78) 7 PP 
Pp r i 


where P_, GG, denote the time-average power input and 
i 


S(r,98,¢) is the time-average power density in the radiated 
field. GPP and GIT denote the power gains associated with 
the phi-polarized and the theta-polarized components of the 


field, respectively. 


The use of the variables JFLAG and KFLAG are described 
in subroutine SGANT. 


CALLED BY: MAIN 

CALLS TO; GDISS 
GFF 
RITE 
SQROT 
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GGMA 


PURPOSE: to calucate the mutual impedance between two 


filamentary monopoles with sinusoidal current distribution. 


METHOD: As stated in subroutine SGANT, the mutual 
impedance of coupled dipoles may be expressed as sum of four 
monopole-monopole impedances. This subroutine calculates the 
mutual impedance with closed-form expressions in terms of 


exponential integrals. 


For skew monopoles it can be shown that the 


monopole~monopole mutual impedance is given by: 


a e e G G 
ij j1 12 22 


where m= 2/1, n = 2/j and 


1) 
16 n sinh d_ sinh d 
1 2 


The functions tre are defined by: 
a 


qz cos § 
7 
F = 2 sinh d e ECR ete0zeeces Oo aargt) 
ik 1 a, al 


k 
where gq = (-1) , 7 and ss are the lengths of the monopoles 


being considered. The functions Ss are defined as follows: 
i 
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G = E(R + qz + g't -—- jqt) + E(R + qz + gq*tt + jaqt 
ie eee ja") + B(R, + az, + g ja") 


re le + q2 fee? te — gee) = aap + 27 g q't + jq") 


i k 
where q = (-1) , q' = (-1) , and q* = gb + q'c with 
b = c cos ¥ and c = d/sin J. The angle 9 is the angle formed 
by the apparent intersection of the two monopoles. This 


will be discussed later in detail. 


ik 
In the. above equation for G  , t denotes the position of 


an observation point somewhere on monopole 2. ae and Be are 


the distances from the endpoints of monopole 1 to this 
observation point. Finally, the E functions are defined as 


follows: 


+iqtt 
a +iq 
where a and q" are real quantities with dimensions of 


length, ais a function of t, = = oe og So = ihe and 


y= jw\/fus. The integral above is evaluated by subroutine 


S2 
T2 


T2 
Se S1 A 


¢ 
¢ 


| 
I 
J 
4’ 


PARALLEL SKEW. 


To explain the input data for GGMM, refer to the above 
figure. If the monopoles are parallel, then the new 
coordinate system is defined such that the new z axis is 
parallel to the monopoles. The coordinate origin may be 


selected arbitrarily. $1 and S2 denote the z coordinates of 


> 1 


the endpoints of the test monopole, T1 and T2 are the 
coordinates of the endpoints of the expansion monopole, and 
D is the perpendicular distance (displacement) between the 
monopoles. The mutual impedance of parallel monopoles is 
calculates in the last part of GGMM below statement 5. 


For skew monopoles, let the test monopole s lie in the 
xy plane and the expansion monopole t in the plane z = D. (D 
is the perpendicular distance between the parallel planes.) 
If the monopoles are viewed along a line of Sight parallel 
with the Z axis, the extended axes of the two monopoles will 
appear to intersect at a point on the xy plane. Let s 
measure the distance along the axis of the test monopole 
with the origin at the apparent intersection. S1 and S2 
denote the s coordinates of the endpoints of the test 
monopole. Similarily, let t measure the distance along the 
axis of the expansion monopole with the origin at the 


apparent intersection. T1T1 and T2 denote the t coordinates 
of the endpoints of the expansion monopole. Let S snd € be 
unit vectors parallel with the positive s and t axes, 
respectively. Then CPSI = Se€ = cos J. The monopole lengths 


are qd andd. 
Ss c 


The output data from GGMM are the impedances P11, P12, 
P21, and P22. In defining these impedances, the reference 
direction is from S1 to S2 for the current on monopole s, 
and from T1 to T2 for the current on monopole t. In the 


impedance P_, the first subscript is 1 or 2 if the test 
1) 


dipole has terminals at S1 or S2 on monopole s. The second 
subscript is 1 or 2 if the expansion dipole has terminals at 


T1 or T2 On monopole t. The monopole lengths d and d.. are 
Ss 


assumed positive in defining the input data CGDS, SGD1 and 


D2 


SGD2. 


CPSI = 1 or 


For parallel monopoles, 
T2 are cartesian coordinates for parallel 
spherical coordinates for skew monopoles. 


monopoles, the radial coordinates S1, S2, 


infinity as the angle ¥ tends to zero or n. 


=, 


a4, 052, T1,-and 


monopoles and 
For skew 
T1, and T2 tend to 


Therefore, if 


the monopoles are within 4.5° of being parallel, they are 


approximated by parallel dipoles. 
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GGS 


PURPOSE: to calculate the mutual impedances between 
two filamentary monopoles with sinusoidal current 
distributions. 


METHOD: The monopole-monopole mutual impedance as 
defined by SGANT is calculated using the equations defined 
in subroutine GNF. The endpoints of the axial test monopole 
S are (XA,YA,ZA) and (XB,YB,ZB), and the endpoints of the 
expanSion monopole t are (X1,¥1,Z21) and (X2,Y2,Z22). DS and 
DT denote the lengths of monopoles s and t, respectively, 
CAS, CBS and CGS are the direction cosines of monopole s, 


and CA, CB and CG are the direction cosines of monopole t. 


The effects of ground for vertical co-linear nonopoles 
are applied ina slightly different manner than mentioned 
previously. As with self impedance calculations, the test 
monopole and the expansion monopole are laterally displaced 
by the wire radius.- This lateral displacement is used to 
determine the angle of incident. This technique is applied 


at statement 8. 


If INT = 0, GGS calls GGMM for the closed form impedance 
calculations. Otherwise GGS calculates the mutual impedance 
via Simpson's-rule integration with the following number of 
Sample points: IP = INT + 1. Lf the monopoles are parallel 
With small displacement, GGS calls GGMM to avoid the 


difficulties of numerical integration. 


Since the point (X,Y,2Z) of subroutine GNF lies on the 
expansion monopole t, T is the integration variable and is 
measured from (X1,Y1,Z21). C11s the current at T for the 
mode with terminals at (X1,Y1,Z1), and C2 is the current at 
T for the mode with terminals at (X2,Y2,Z22). C denotes the 


Simpson's-rule weighting coefficient. 
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Below statement 7, GGS performs some analytic geometry 
in preparation for calling GGMM. The remainder of this 


section is concerned with this preparation. 


Let s denote a unit vector in the direction from 
(XA,YA,ZA) toward (XB,YB,ZB). Also let € denote a unit 


vector from (X1,Y¥1,21) toward (X2,Y2,2Z2). Then Sef = 
cos 8 = CC where 98 is the angle formed by the axes of the 


two monopoles. Let monopole s lie in one plane P and 
S 
monopole t lie in another parallel plane ae CAD, CBD and 


CGD are the direction cosines of the unit vector d = € x5 / 
Sin @ which is perpendicular to both planes. To obtain the 


distance DK between the two planes, a vector Rk 1s 
constructed from (XA,YA,ZA) to (X1,Y1,Z1) and take DK = 


Red. 


ii 


A line is constructed from (X1,Y1,2Z1) to the test 
monopole, such that the line is perpendicular to the test 
monopole. SZ denotes the s coordinate of the intersection of 
this line with the test monopole, and the cartesian 


coordinates of this intersection are XZ, YZ, and ZZ. The 


direction cosines of s x d are CAP, CBP, and CGP. 


From the point (X1,Y1,Z21) in plane ff a line is 


constructed perpendicular to the point (XP1,YP1,ZP1) in the 
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llel with d and has length DK. 


1S para 


- This line 


plane P 


Let R represent a vector from (XZ,YZ,2Z) to (XP1,YP1,2P1). 


broutine 


in su 


S71 and T1 are defined 


P1 denotes R (S x dq). 


GGMM. 


SGANT 


CALLED BY 


GGMM 
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PURPOSE: to calculate the near-zone electric field of 


a Sinusoidal electric monopole. 
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METHOD: An electric line source is located on the z 


axis with endpoints at Z and z, as Shown in the above 


figure. Let the electric monopole have the following current 


distribution: 


I sinh Y(d - 1) + I sinh yl 
= (1) = a | 2 


Sinh yd 
where I, and I, are the endpoint currents, Y 1s the complex 


propagation constant of the mediun, d = z, 7 2, 1s the 


source length. The cylindrical components of the field are 


E(8) = 0 and 


n -YR “YR ; 
E(@€) = C (1, = 1-1 ee 2) sinh yd 
Yn€@ sinh Yd 2 
-YR 
+ (I, cosh yd - I) e 1! cos oT 


-yYR 
t (I. cosh yd - I) e / ¢ cos coe 
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= R 
E(z) = ae ee f({I - I cosh Yd) 3 f 


4n sinh yd 1 Z 
+ I Hadj een 
(I. . cosh yd) e . ] 
1 
where is the intrinsic impedance of the medium and where 


(©,@,z) denote the cylindrical coordinates in a coordinate 


system centered at the endpoint of Zo: 


These expressions exclude the field contributions fron 
the point changes at the endpoints of the line source, since 
these charges disappear when two monopoles are connected to 


form a dipole. 


Let the coordinate s measure distance along the test 
monopole with the origin at (XA,YA,ZA). From any point 
x5 ¥e2, a line is constructed perpendicular to the monopole. 
SZ denotes the s coordinate of the intersection of this line 


with the monopole. The length of the line is the radial 


2 
coordinate @, and RS denote @ . R1 and R2 are the distances 


from (XA,YA,ZA) and (XB,YB,ZB) to the point (X,Y,2Z). 

In the statements above statement 1, the above equations 
are solved; and after statement 1, the cartesian components 
(E ,E ,E ) of the field are determined. If a ground plane is 

a“ y Zz 


present (IGRD>0O) the reflection coefficients are applied to 


the cartesian components before returning to the calling 


program. 
CALLED BY: GNFLD 
CALLS TO; NOWE 
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GNFLD 


PUKPOSE: 


intensity at a given point. 


METHOD: 


field of each wire segment, 


obtain the near-zone field of the wire antenna. 


in a manner similiar to FI of subroutine SGANT. 


loop currents calculated by 


to calculate the near-zone electric field 


This subroutine calls GNF for the near-zone 


and sums over all segments to 


FI is used 


CJ(I) is the 


subroutine GANT1. 


The use of the variables JFLAG and KFLAG are described 


in subroutine SGANT. 


CALLED BY: MAIN 
CALLS TO: GFF 
L 
SUBROUTINE GNFLO (IAs By sENM, DLP LZE 19 2MDe Ns ~ND, NM, AMeCGDy SG0, ETA,GAM L 
PyEX 
eevee sey tL Peso evi. Fee Ree ez ETALGAM 3 
COMPLEX ERK : 
SEU Eom atl Poa R eae OF a, B2CRDe R3BCDD» OCLDy MCLs VERDe 201 6 
1) 7 
MENSION MOCINM, 4), NOGI) : 
DATA PA 8 tents b. 283187 d 
Ex) = (5090) 10 
EY = {-Q,-0) 12 
El = (.0,.0) 3 
: OO 2 Kz{,NA is 
KA = TACK) ts 
KA = IB(K) 3 
Age er IGRO=-1 62 eE6D ‘8 
WA, ° Al DeXPe VP, lPeAMe DI Ke 
F Kade VYEKAD?LIKA (KBE eVEKBE, 2EKB 0 
Br ea eee EAM EXL EY LeELIeEXSeEV2rEL2e IGRO,ERRS 20 
GRO = NGRO 21 
NOK = NO(K) 33 
C 
00 2 11#=(_,NOK 24 
= mocks Li) 36 
DE (KB.EQ.12499) GOTO 1 27 
NF (KB. EQ. INC ID) Fle-l. 28 
Exe EXeFIeEXleCscd) 29 
Ey e EYOFLOEYIECIJOD) 30 
El = ELOFLSEZLECIOI 3: 
GO TO 2 Z 
1 TE (KAsEQ.I340DD Ebene 33 
EX = EX@FLPEXZ°C ICL) 3% 
EY = EVeFIsEYZeCI{ 1} 39 
El = ELOFIOEL2Z°C SE 36 
2 CONTINUE 38 
C 
RETURN 23 
ENO 
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LEFT 


PURPOSE: to determine position {location) of the Left 
paren symbol on the input data card. 


METHOD: The character search begins in the column 
passed to the subroutine. On returning to the calling 
program the argument passed is the column following the left 


paren symbol. 


CALLED BY: READ 


CALLS TO; NONE 


eRe 2 OD Ow We 


= es 
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in 
In 


is inserted 3 


haracters on the polar plot. 
the grid numbers labels 


1a Cc 


the statements above statement 4. 


— Or 


in 


inser 


The peroid character ({ISM(2)) 
the proper position 


POLPLOT 


to 
NOWE 


the statements after statement 4, 
are inserted on the horizontal axis. 


LINECK 
PURPOSE 
METHOD 
CALLED BY 
CALLS TO 
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NUMB 


to place degree numbers on the polar plot. 


PURPOSE 


The current line which 1S being printed is 


McoTHOD 


passed to the subroutine in the calling argument. 


If this 


these numbers are placed in 


line contains degree numbers, 


the correct position by the IF statements. 


PTPLOT 


CALLED BY 


NONE 


CALLS £0 


ANIA PAO OPO MN AFA OM DFO MU AEN O™ DROMRNM ENO DAO ANMEN 


SUBRUUT ENE NUMB (Y} 


THES SUBROUTINE PUTS DEGREE NUMBERS ON POLAR GRID 
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NUMBER 


PURPOSE: to convert alpha-numeric numbers to floating 


or fixed point numbers. 


METHOD: After initially determining the sign of the 
number, the DO LOOP ending at statement 6 Scans each 
character beginning at N1. The DO LOOP ending at statement 
3 terminates the outer DO LOOP if the character being 
compared is not an alpha-numeric number. The DO LOOP ending 
at statement 5 converts the alpha-numeric number to an 
actual number. Below statement 7, the multiplier correction 
is applied to the floating point number before returning to 


the calling progran. 


CALLED BY: READ 
CALLS TO: NONE 
SUBROUTINE NUMBER INI gN2eXe 1X) l 
COMMUN /A/ ALB0) 2 
CATANGTMORmREL | wale 3 eae +5" : | e ye ¢ge toes 3 
Data AMNUS PLUS, POINTsI-*, tee, teyy 5 
chien ee oo Bead 6 
NSIGN © 0 4 
i= -] 9 
ix 20 10 
1SET © D 11 
LF [TAQN).EQePLUS) N=NOL 12 
IF [AIN).NE.AMNUS) GO TO 1 . 13 
NSIGN = 14 
: N © Nel 15 
1 00 6 L=N,80 i$ 
[F TACL)CNESPOLNT) GO To 2 18 
ISer = | 19 
GO TC 6 20 
co? FR OUTSEr.€.10 11 © 1161 21 
00 3 K=1el0 $e 
IF (A1L)+EQ.8(K)) GO TO 4 54 
¢ 3 CONTINUE 35 
Go to 7 36 
C 56 
& DO 5 K=1,10 29 
KK © K-1 30 
IF 1ALL}eEQeB(K)) NUMBEKK 31 
co CONTINUE 32 
1X © NUMB61001Xx 3a 
N2 = le] 35 
< & CONTINUE 36 
? 
TAF INSIGN.EQ.1) 1K = -1X 36 
IF E11et1T.0) Ll © D 20 
X = Y/L 108911) 41 
1F LA(N2) -EQ-POLNT) N2"N2¢1 4&2 
If {AINZ).EQ.aK} X = X#1000, 43 
le (AUN2)-EQ.aN} X © x20.00f 44 
IF (A(N2}.EQ.AU) X © Xe0.000001 4&5 
LE CCAINZ) EQ. AK) SOR. (A(NZ) «EQ+AM) ORTAINZ) .EQ2AUI) NZONZ¢L 46 
RETURN <5 
ENO 49 


67 


POLPRT 
PURPOSE: to control the plotting of the polar plot. 


METHOD: This subroutine is the main subroutine in the 
polar plot package and is responsible for calling the 


various subroutines of the package. 


The scale factor, S, must be changed according to the 
printer characteristics. The scale factor in this subroutine 
is set for ten, 10, characters per inch for the abscissa and 
eight, 8, characters per inch for the ordinate axis. 
Therefore S = 10./8. 


After initializing DATAX, DATAY, and X, the input data, 
Y, is scanned to determine the normalizing factor. If this 
normalizing factor is less than 1.E-32, an error statement 


is printed and the plotting is abortted. 


In the DO LOOP ending with statement 8, each line of 
the polar plot is printed after a call is made to PTPLOT to 
establish the ploar grid information. The variable, DIM, is 
used to as a scaling factor for the polar plot. The value of 
1.0 will cause all of input data to be plotted, however, if 
only the values less than one-half of the normalizing factor 
are of interest, then DIM can be set to .5. This will 


enlarge of the center of the polar plot. 


CALLED BY: MAIN 
CALLS TO: PTPLOT 
SART 
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PTPLOT 


PURPOSE: to establish the grid information for the 
polar plot. 


METHOD: In the DO LOOP ending at statement 1 the 
alpha-numeric characters are transferred to ISYN in order 
to pass via COMMON to other subroutines. In the statements 
following statement 2, the equations for the plotted 
concentric circles are established. Below statement 7 the 


grid marks on the 090-270 axis are inserted. 


CALLED BY: POLPRT 
CALLS TO: LINECK 
NUMB 
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READ 


PURPOSE: to interpret and translate the input data 
cards. 
METHOD: The program utilizes free format for the data 


cards, that is, the program uses character recognition to 
determine which parameters are being read. In the IF 
statements containing A(1), A(2), A(3), and A(4), the first 
four characters on the data card are compared to the first 
four letters of the key words. This will determine the type 
of parameters that card contains. The other IF statements 
determine which parameters are being read. 

Subroutine BLNK is called to remove the blank spaces on 
the parameter cards. Subroutines EQUAL and LEFT are called 
to determine the position of the equal character and the 
left paren, respectively, Subroutine NUMBER is called to 
convert the alpha-numeric characters to numbers, either 
fixed or floating point. This numerical value is assigned 


to the parameter just determined. 


A detailed explanation of the data cards is found in 
appendix II titled "USERS MANUAL". 
CALLED BY; BLNK 
EQUAL 
LEFT 


NUMBER 
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RITE 


PURPOSE: to generate a list of branch currents fron 


the input loop currents. 


METHOD: The generation of branch currents is 
accomplished in the DO LOOP ending at statement 2. The 
branch currents are stored in CJ(I) by the latter part of 
the DO LOOP ending at statement 3. If the branch currents 
are reguested for output (IWR positive), the DO LOOP ending 


at statement 5 accomplishes this. 


CALLED BY: GANT1 
. GFFLD 
CREES TO: NONE 


80 


NPA OM DP OAM FAO PO NMEA OM DAO MANNE AO DHOMNM ES AOD 
mt omel Pt ome ond od (CIEE AGG AI ENE NE NEE 9 PP PP NCA OD bo oe 
-_— vam omel Pre) stoma ™ PITS P ES TOAAMAMMAMAMAMAA OOO OOOO 


mm 
~~ xX 
o = 
& sos 
— uw ~ 
— = o > 4 
oOo = =x «6S 
ae © = eo 
& td 
2 = ao Wi ed 
. om am . « 
oO = am & — ww 
WJ <o™ ~~ « 
to od =] IW Ww © 
i mee WW o ~ 
-— = a ad —< Po 
z °& Ld Ld AD Wee 
- x aeaw eo Wu 
© eee =~ Ou 
ze OD ~ O« 
-_ — <o™ w Fr 
oOo «= =) “HH o a 
, WIL) o am f& @ 
e * — <7 —< — & 
mo Oox - NM 
aed <<q— “Ww ery 
= bad > es) * *® @ 
Nn = omg om ed —“ Me 
— q<ae me LL 
~ = ome nee Ww WO 
~—=a NN NN a 5 
— om ZzZe <x! he 
ba i eae — wo 
ace : bo mm) - -—Zz 
z= <aq< “w wo 
ay eae . Wa 
*D— oog = ox 
xr“ pra he 3 — © @ 
zw 559 N IAL) w Oe 
—. = Ww . PAW - W»« 
oie ° ° <n ° x NN & ay 
ol Qo lm oe ~ 2 o od N we 
ons ms ' mC 2. w hh = we 
pgeamead a s << <f =< gq AWA NM @w 
a—~-mD © = = xx OQ « e —= FO 
os oma Pat OoOu~ um a< -_ xX rh) soo e* a> 
othe = oy < x x =x Qe 
Wi) am mm ame ate el E Oo - < ~ anand —w or 
— od wma) et) dia < ° ~ 5 D « = N “5 
we —_ 4 ww we =A 7ZxX ~ © oO U I< an am & e «<a 
<< oo OF Nee Me OvLa4 oO -« oo Oe VEX ° ox ~N OX 
= x ce <= Pl © eel << “ww © © wus —- Fm yv'i — ig oo —— fe 
we Zon eo ox = a oe > Wueewr - gD OO) Fx wo vw oom = Qe we 
ZoZzo *x x OOw =—=es COO Oc WE DovwH 9 oe e ome Xo ff oda Do b= be 2) xx ON 
- Ge maw 8 0 tw ews wow Lida o ow 7 MOTO =UZFZ0UVU i775 WW eo e ay Ve | 
m— fe ome WG CX) nue woe —w © © oN eWD 1? WWD eg de ra HuUeU WWW e eH & ww ww & 
Dwye 2 on ome “—O=Dn ft 2 ~~ WZ zzx 8 xx 8 WYSOII~ ~ - 
Cuz nw X xxNOxwue ~ wendoe mo et . 8 « ¢ OU z= -- => 
eoux mW w NN haw wwe te wux wer sr be ~ 3s w= Yo wou <ow eT. 4 <q <a 
ertad <aox a < oF ~~ ~9>5 2 - =o «05 Oo uvrewe r-D CX Xe 
Boas a io BS = “au Ou 3o WUeOvUvUUUds uWwudgte Onx 70 Coo wat = x w oO S22 
mk ret OL WO <f Ot mC oe ms bat ee 
== = z= Sr liad: Wt OOO —- E za Le LL om. CD 
—= ON Laat] i J — od 
a om @ 


34 


SAL 
PURPOSE: tomsobtedata for polar pet: 


METHOD: This subroutine sorts the values of the 
points to be plotted by the polar plot package starting with 
the greatest positive value of y to the greatest negative 
value. In the DO LOOP ending at statement 1, the value of 


(x -y.) 1S interchanged with the value of (x ,y ) if y. is 
el Hs 4, 
greater than y.. 
1 


CATED BYs BOLE RE 


CAEES TO: NONE 


SUBROUTINE SART (OATAX,OATAY,N) " 
; DIMENSION DA&TAX(500), OATAY (500) 2 
c THIS ROUTINE SORTS OATA IN OATAY BY MAGNITUOE r 
; NN = N= ; 
OO 2 I=1yNN ; 
NM © Jel 9 
C 10 
O00 | J=NM,N 11 
IF (OATAYLI).GE.OATAY(J)) GO TO 1 12 
STOR = DATAY(L) 13 

DATA YCT} = DATAY(YJ) : ie - 
DATA YCJ} = STOR 15 
STOR = DATAXCL) 16 
OATA XC1} © QATAX(J) 17 
OATA X(J} = STOR 18 
1 CONTINUE 19 
2 CONTINUE s? 
RETURN $$ 
ENO 24 
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SGANT 


PURPOSE: to calculate the mutual impedance between 
filamentary monopoles. 


2 Z 
32 coeceee 2S arc : 
P12, ~* 
ea ~7” 
1 Swe 1 
o* Se 2 
P21° 5 
oe" Say 
BO | eweeeee== P2 2ee=2===s= D 
1 1 
SOURCE RECEIVING 
DIPOLE DIPOLE 


Z= P11 + P12 + P21 + B22 


METHOD: In the induced emf formulation, the nutual 


impedance of coupled dipoles is 


Z=- i (t) E(t) dt 


where Bw) denotes the current distribution (normalized to 
unit terminal current) on dipole 2, and ae is the field 


of dipole 1 when it transmits with unit terminal current. 
Distance along the axis of dipole 2 is denoted by the 


coordinate t. a May be expressed as the sum of the fields 


from each of the monopoles comprising dipole 1. Furthermore, 
the integral is the sum of the integrations over each of the 
monopoles comprising dipole 2. Thus, the dipole-dipole 
mutual impedance may be expressed as the sum of four 


monopole-monopole impedances. 


It may be convenient to draw the above figure in terms 


of monopoles with the current distribution shown as dotted 
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lines. (The monopole letters remain the same.) 


The surface impedance 1S calculated just above statement 
2. BO1 denotes aye / a where ur and a are the Bessel 


functions of order zero and one with complex argument, ZARG. 
It is assumed that all the wire segments have the same 


radius, conductivity and surface impedance. 


In the DO LOOP ending with statement 3, SGANT calculates 
the segment lengths D(J). DMIN and DMAX denote the lengths 
of the shortest and longest segments. If the wire radius or 
the segment lengths are clearly beyond the range of 
thin-wire theory, N is set to zero at Statement 4 followed 


by RETURN to the main program to abort the calculation. 


At statement 5, the program selectS a Segment K, anda 
few statements below this it selects another Segment L. K is 
a segment of test diple I, and L is a segment of expansion 
mode J. The mutual impedance between segments K and L is 
obtained by calling subroutine GGS or GGMM. In statement 18, 
this impedance is lumped into C(MMM). The mutual impedance 
a5 between dipoles I and J is the sum of four 


segment-~segment impedances. 
The variables IFLAG and JFLAG are used if a ground 


plane 1s present for the calculation of the mutual impedance 
elements. If IFLAG 1s egual to JFLAG, the mutual impedance 
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terms will not have the effects of a ground plane since both 
monopoles lie on the same side of the ground interface. If 
the monopoles are on the opposite sides of the interface 
(IFLAG not equal to JFLAG), the reflection coefficient 
correction must be applied to the mutual impedance elements. 
This same technique is applied in subroutines GNFLD and 
GFFLD. 


In SGANT, segment K has endpoints KA and KB, and segment 
L has endpoints LA and LB. It is convenient to think of KA 
and KB as points 1 and 2 on segment K, and LA and LB as 
points 1 and 2 on L. The four segment-segment impedances 
can be defined as P(IS,JS). The first subscript IS refers 
to the terminal point on segment K, and the second subscript 
JS refers to the terminal point on L. Thus IS=1 or 2 if 
dipole I has its terminal point I2(I) at KA (point 1) or KB 
(point 2), respectively. Similarly, JS=1 or 2 if mode J has 
its terminal point I2(J) at LA or LB. The impedances 
P(IS,JS) are defined with the following reference directions 
for current flow: from point 1 toward point 2 on each 
segment. If dipole I has this same reference direction on 
segment K, FI=1; otherwise FI=-1. Similarly FJ=1 or -1 in 
accordance with the reference direction for mode J on 
segment L. In statement 18, P(IS,JS) is multiplied by FI and 


FJ before its contribution is added to Z . 
13 


Subroutine GGMM calculates the impedances Q(KK,LL) which 
are like the P(IS,JS) but have different conventions for 
reference directions and subscript meaning. The 
transformation from the Q impedances to the P impedances is 


accomplished in the DO LOOP ending with statement 13. 
If the wire has finite conductivity, the appropriate 


modification is applied to the impedance matrix just above 


statement 15. The terms arising from the dielectric shell 
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on an insulated segment are obtained from subroutine DSHELL 
just above statement 16. Finally, the lumped loads, ZLD, are 
added to the diagonal elements of the impedance matrix in 
the DO LOOP ending at statement 23. 


K is a segment of test dipole I, and L is a segment of 
expansion mode J. When the segment numbers K and L are 
equal, SGANT calls GGMM to obtain the mutual impedance 
between two filamentary electric monopoles. These monopoles 
are parallel and have the same length. Monopole K is 
positioned on the axis of the wire segment, and monopole L 
is on the surface of the same wire segment. Thus, the 
displacement is equal to the wire radius. The two monopoles 


are Side~by~side with no stagger. 


When segments K and L intersect, SGANT again calls GGMM 
for the mutual impedance between the two filamentary 
monopoles. Monopole K is situated on the axis of wire 
segment K, and monopole L is on the surface of wire segment 
L. The axes of segments K and L define a plane P,and 
monopole K lies in this plane. Monopole L is parallel with 
plane P and is displaced from it by a distance egual to the 


Wire radius. 


CALLED BY: MAIN 

CALLS TO: CBES 
DSHELL 
GGMM 
GGS 
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SORT 
PURPOSE: to define 
METHOD: In the DO 


of dipoles is defined by 
(the endpoints of dipole 


dipole I); and the vectors JA(I) and 


the set of dipole modes. 


LOOP ending 


filling the 


I); 


I2 (I) 


( 


at statement 3, the set 
vectors I1(I) and I3(I) 
the terminal point of 
JB(I) (the monopoles 


comprising dipole I) with the node numbers and segment 


numbers, respectively. The DO LOOP ending at statement 8 


determines MD(J,K) (the list of dipoles sharing segment J) 


and ND(K) (the number of dipoles sharing segment J). 


CALLED BY: MAIN 


CALLS TO; NONE 
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SQROT 


PURPOSE: to solve the set of simultaneous equations to 


determine the currents on the thin wire structure. 


METHOD: This subroutine considers the matrix equation 
ZI = V which represents a system of simultaneous linear 
equations. NEQ denotes the number of simultaneous equations 


and the size of the matrix Z. 


On entry to SQROT, S is the excitation column VY. On 
exit, the solution I is stored in S. Z({I,J) denotes the 
symmetric square matrix. Also on entry, the upper-right 


triangular position of Z(I,J) is stored by rows in C(K) with 


K = (I - 1) * NEQ - (I * I) S2+9. 


If 112 = 1, SQROT will transform the symmetric matrix 
into the auxiliary matrix (implicit inverse), store the 
result in C(K) and use the auxiliary matrix to solve the 
Simultaneous equations. If I12 = 2, this indicates that C(K) 


already contains the auxiliary matrix. 

The transfromation from the symmetric matrix to the 
auxiliary matrix is accomplished in the DO LOOP ending at 
statement 5. The solution of the simultaneous equations is 


accomplished in the remainder of the progran. 


CALLED BY: GFFLD 


CALLS TO: NONE 
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SYMBOL DICTIONARY 


characters of the input data cards 
backscattering phi plane angle for plotting 
backscattering theta plane angle for plotting 
bistatic scattering phi plane angle for 


plotting 


bistatic scattering theta plane angle for 
plotting 


absorption cross section for phi polarization 


absorption cross section for theta 
polarization 


far-zone phi plane angle for plotting 
far-zone theta plane angle for plotting 
radius of the thin wire of the structure 


outer radius of the dielectric shell of the 
insulation of the wire 


elements of the open-circuit impedance matrix 
branch currents for the structure 

cosh yd for a given segment 

loop currents for the structure 

conductivity of the wire 

length of a given segment 

extinction cross section for phi polarization 


extinction cross section for theta 
polarization 


Cadiation efficiency 

loop currents induced by a phi polarized wave 
phi-polarized far-zone feild of the dipole 
mode 

scattered electric field in the phi direction 
due to a phil polarized wave 


scattered electric field in the theta 
direction due to a phi polarized wave 


complex permittivity of insulation 
complex permittivity of ambient medium 
complex permittivity of ground 

EP4/EP3 


relative dielectric constant of insulation 
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ER3 


ERY 
ET 


ETA 
ETE S 


BIT 


ETTS 


EX 
EY 
EZ 
E0 
FHZ 
FMC 
GAM 


GG 
GPP 


GTT 


HGT 
IA 

iB 

1E)e ES © 
ICARD 
ICJ 


IFLAG 
IGAIN 
IGRD 
INC 


INEAR 
INM 


INT 


relative dielectric constant of the ambient 
medium 


relative dielectric constant of the ground 


loop current induced by a theta polarized 
wave 


intrinic impedance of ambient nediun 


scattered electric field in the phi direction 
due to a theta polarized wave 


Cae polarized far-zone field of the dipole 
node 


scattered electric field on the theta 
direction due to theta polarized wave 


near-zone electric field in x direction 
near-zone electric field in the y direction 
near-zone electric field in z direction 
8.854E-12 

frequency in hertz 

frequency in megahertz 


intrinic progration constant of the ambient 
medium 


time-average power input 


power gain associated with the phi polarized 
component 


power gain associated with the theta 
polarized component 


height of the structure above ground plane 
first node of a given segment 

second node of a given segment 

indicator for bistatic scatter calculations 
indicator for the data cards 


dimension corresponding to the number of 
Simultaneous linear equations 


indicator for program termination 
indicator for antenna gain calculations 
indicator for presence of the ground plane 


indicator for the type of far-zone 
calculations 


indicator for near-zone calculations 


dimension corresponding to the number of 
monpoles 


number of integration steps 
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ISc 
ISCAT 
IWR 


I12 


indicator for 
indicator for 


indicator for 


the insualition 
backscatter calculations 


current distribution output 


endpoint node of a given dipole 


indicator for 
terminal node 


auxiliary matrix 


number of a given dipole 


endpoint node number of a given dipole 


first segment 


number of a given dipoile 


second segment number of a given dipole 


print indicator 


list of generator/feed locations 


indicator for structure load 


list of impedance/load locations 


maximum of the number of segments connected 
to any one given node 


list of dipoles Sharing a given segment 


minimum of the number of segments that 
connected to any one given node 


indicator for error printout 


number of simultaneous linear equations 


total number of dipoles sharing a given 
segment 


indicator for antenna calculations 


number of segments 


indicator for polar plot 


angular frequency 


phi angle for far-zone calculations 


scattering cross section for phi polarization 


scattering cross section for theta 
polarization 


Sinh Yd of a 
conductivity 
conductivity 


conductivity 


echo 
echo 
echo 


echo 


area 


area 


area 


area 


given segment 

of insulation 

of the ambient medium 

of ground 

phi incident-phi scattered wave 
phi incident-theta scattered wave 
theta incident-phi scattered wave 


theta incident-theta scattered wave 


oF 


TD2 
TD3 


loss tangent of the insulation 

loss tangent of ambient mediun 

theta angle for far-zone calculations 
2n (6.28318) 

1.2566E-6 

antenna complex driving voltages 

list of VG's 

x-coordinate of each node 

list of XP's 

x-coordinate for near-zone calculations 
y-coordinate of each node 

list of YP's 

y-coordinate for near-zone calculations 
complex power input 

z-coordinate of each node 

complex load at a given node 

list of ZLD's 

list of ZP's 

z-coordinate for near-zone calculations 
surface impedance of the wire 


antenna input impedance 
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USER'S MANUAL 


The Antennas~-Scatterers Analysis Program (ASAP) for thin 
wire structures in a homogenous conducting medium performs a 
frequency domain analysis of antennas and scatters. The 
program is applicable in the presence of either a perfect or 
a finite ground. This appendix will describe and explain 
the data cards necessary to execute the compute progran. 
Although the program was written for the IBM 360 computer 
system, it can be executed on another system with minor 


modifications. 


The program utilizes piecewise sinusoidal expansion for 
the current distribution with Kirchhoff Current Law enforced 
everywhere on the structure. If the structure contains end 


points, the currents at these points are assumed to vanish. 
I. Program Limits 


The thin wire assumptions are questionable and the 
accuracy and convergence deteriorate if the radius of wire 
utilized for the structure exceeds 0.01 of a wavelength, if 
the longest segment is greater than one~fourth of a 
wavelength, if the length ratio of the longest and shortest 
segments exceeds 100, or if the total wire length is less 
than 30 times the wire diameter. If a wire is bent sharply 
to form a small acute angle (less than 30 degrees), the thin 
wire model is guestionable. It 1S assumed that the wire 
conductivity greatly exceeds the conductivity of the ambient 
medium. For insulated wires, the dielectric layer is assumed 


to be electrically thin. 
II. Minimum Data 


The minimum data necessary to execute the program is: 
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a. description of structure 


b. radius of wire used for the structure 


The program will default to the other parameters necessary. 


The default parameters are: 


a. wire for the structure is copper 


b. frequency of operation is 300 mhz 


c. homogneous medium is free space 


A more detailed explanation of the defaults will be 


discussed when the data card for the parameter is described. 


IIIT. Outputs 


In antenna problems, the output includes structure 
currents, impedance(s) of feed(sS), gain, polar radiation 
plots, and near field calculations. In bistatic scattering 
problems, the output includes structure currents, complex 
elements of the polarization scattering matrix, polar 
reradiation pattern plots, and echo areas produced by a 
plane wave. For backscattering problems the output includes 
absorption, scattering and extinction cross sections in 
addition to the outputs of bistatic scattering. Most of the 
Outputs are Suppressed and must be requested. Since the 
program can produce a large volume of output, care should be 


exercised until the user is familiar with the outputs. 
IV. Data Cards 

The Analysis Program utilizes free format for the data 
cards, that is, the program utilizes character recognition 


to determine which parameters are being read. Data placement 


(location) on the input card is not critial. Blank 
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characters, on all input cards but the COMMENT data card, 
are ignored and may be used at the discretion of the uSer. 
Since character recognition is used, only the first four 


characters of the key words must be present and correct. 


The format for the COMMENT CARD utilizes standard 
FORTRAN format (i.e. 'C* in column 1 followed by at least 
four blanks). The COMMENT CARD is the only type of input 
card that position in the data block is critial. This 
(these) card(s) must be placed at the beginning of a data 
block. A data block is a series of related data cards. 
Several data blocks may be used to define an analysis 
problem. This will become clear when the termination cards 
(END, STOP, or CHANGE) are discussed. There is no limit to 
the number of comment cards that may be used. AS a check for 
the user, all input data cards will appear on the output as 


they appear in the input deck. 


The format of other data can be of one of two forms: 


ae type of card (option 1/foption 2/....-..--) 


b. parameter (value) . 


The type of format to use will be apparent as the individual 


data cards are discussed. 


The numerical values for the parameters may be stated in 
any one of the following forms. The program will translate 
the number to the proper form for the specified parameter, 
either fixed or floating point. All of the following 


examples have the same value. 
0.0001 or .0001 or 100.U or 100U or .1M or O.1M or .0000001K 


=6 = 3 
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1. WIRE This card is used to define the parameters 
associated with the wire utilized by the thin wire 


structure. Two options are available and are defined as: 
RADIUS=value of the radius of the wire in meters 
CONDUCTIVITY=value in megamhos per meters . 


The wire data card must appear in the first data block to 
define wire radius. The default value of the conductivity 


is 50 megamhos/meter (copper). 





IRE( RADIUS=.001/ CONDU=28.5) 
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2. INSULATION This card is utilized to define the parameters 
associated with the insulation of the wire used for the 
structure to be analyzed. If this card is omitted, the 
program assumes that the structure is uninsulated. Four 


options are available and are defined as: 
RADIUS=value of outer radius in meters 
CONDUCTIVITY=value in micromhos per meter 
DIELECTRIC=value of relative dielectric constant 
LOSS TANGENT=value . 

The conductivity and either the relative dielectric constant 


or the loss tangent (but not all three) options may be 
stated. 





INSULATION ( RADIUS=.015/ COND=7./DIEL=5) 
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3. EXTERIOR MEDIUM This card is utilized to describe the 
homogeneous medium surrounding the structure. If the mediun 
is free space, this card may be omitted. Three options are 


available and are defined as: 
DIELECTRIC=value of relative dielectric constant 
CONDUCTIVITY=value in micromhos per meter 


LOSS TANGENT=value . 


AS with INSULATION card state either conductivity or loss 
tangent. 


XTE(LOSS=.45) 
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4. DESCRIPTION This card is utilized to describe the shape 
of the wire structure to the program. The user must divide 
the wire structure into segments of the appropriate length 
and number each node starting at one. A node is a point 
where a segment begins or ends. A maximium of four segments. 
can meet at any given node. An isolated wire must contain 
at least two segments and three nodes. The structure is 
described by stating the node numbers that each segment 


connects. The description of a square loop might appear as: 
DESCRIPTION (1-2/2-3/3-4/4-1) . 

The description of a dipole and reflector might appear as: 
DESCRIPTION (1-2/2-3/3-4/4-5 /6-7/7-8/8-9/9-10) . 


If the description will not fit on one data card continue on 
the next card as if the previous card were longer. The 


dipole example might appear as: 
DESCRIPTION (172/2-3/3-4/4-5/ 


6-7/7-8/8-9/9-10) . 


Note that the last character on the card to be continued is 
a Slant (/). As many cards as necessary may be used. The 
maximum number of nodes permitted is fifty. If ground plane 
is present, the maximum number is twenty-five. If a ground 
plane is present and the structure touches the ground plane, 
the lowest node numbers MUST be used for the touching nodes. 
That is, if the structure touches the ground plane at two 
points, node numbers 1 and 2 MUST be assigned to these 


nodes. 





DESCR (1-2/2-3/3-4/4-1) 
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5. GEOMETRY This card is used to state the physicalslocation 
in rectangular coordinates of each node of the DESCRIPTION 

CARD . The rectangular grid is in units of meters. If node 1 
is located at x1,y1,z1 and node 2 at x2,y2,z2 and node 3 at 

X3,y3,z3,etc., the GEOMETRY CARD might appear as: 


GEOMETRY (x1, ¥1-Z21/X2, 12, Z2/X399Y 3723/7. oe oe ee 


As with the DESCRIPTION CARD, continuation cards are 


permitted. 





GEOM(.1,0,-1/-.1,0.1/-. 1,0-.1/.1,0,-.1) 
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6. FEED For antenna analysis the feed point(s) and 
voltage(s) must be stated. In the foremention dipole and 
reflector example if the feeds were at node 2 with a voltage 
source of .5 at an angle of -90 degrees and at node 4 with a 
voltage source of .5 at an angle of +90 degrees the FEED 
CARD might appear as: 


FREDE 2g. 2-59 0748 ,.5,*90) 


The order of the information for each voltage source is node 
number, magnitude, and phase angle. This order is repeated 
until all sources are stated. If the source information will 
not fit on one card, use another card similiar to the 
initial one; that is, repeat the word "FEED". If only one 
voltage source is applied to the structure, only the node 
number must be stated. In the dipole example, if the drive 


is at node 3, the FEED CARD might appear as: 
FEED (3) 


A default source of one voltage at zero degree phase is 
assumed. Voltage sources should only be stated for nodes 


with only two segments. 





FEED(2,-5,-90/4, .5,+90) 
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7. LOAD This card is used to describe the loads to be placed 
at various locations on the structure. The format for this 
card is Similiar to that of the FEED CARD, that is, the 
word "LOAD" is used in the place of "FEED". The order of 
the information on the card is the same. Since this card is 
frequency dependent, it must be changed if the frequency of 
Operation is changed. No default parameters are available. 
The structure is assumed unloaded unless this card is used. 
Once the structure is loaded, it will remain loaded for the 
remainder of the data block series. To unload the structure 


the following card may be used: 


LOAD (-1) 





LOAD (1,120,-45/3,120, +45) 
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8. OUTPUT This card is used to request output data. Most of 
the output is in tabular form. More than one OUTPUT CARD 
is permitted per data block, but not for the same type of 
output. If only the antenna input impedance, antenna 
efficiency, or time-average power input is of interest, no 
OUTPUT CARD is necessary. These parameters are automatically 
printed if a FEED CARD or GENERATOR CARD is utilized. one 
or more of the following options may be used to request the 


various outputs available. 


FAR FIELD=phi initial, phi final, theta initial, theta 


final 


This option gives the components of the electric field 
intensity in the far field as phi and theta varies between 


limits specified in one degree divisions. 


BACKSCATERING=phi initial, phi final, theta initial, 
theta final 


This option gives the absorption, scattering, and extinction 
cross sections, and the complex elements of the polarization 
Scattering matrix for an incident plane wave illuminating 
the structure from the spherical direction of phi, theta as 


both vary between limits specified in one degree divisions. 


BISTATIC=phi inital, phi final, theta initial, theta 


final 


This option gives echo area and the complex elements of the 
polarization scattering matrix for an incident plane wave 
illuminating the structure from the spherical direction phi, 
theta final of the backscattering output option, reradiated 
in the phi, theta direction as both vary between limits 
specified in one degree divisions. A bistatic output request 
must be accompanied with a backscattering request in the 
same data block. 


STEP=value in degrees 


This option will cause any of the above output options to be 
stepped at a different interval size. That is, if one of the 
above options is to be stepped at ten degrees intervals, use 


this option. This option overrides the one degree stepping. 
CURRENT 


This option gives the currents on the structure which are 
produced by the feed/generator voltages and/or the incident 


plane wave of the backscattering request. 
NEAR=x1,y1,2Z1 
Or 
NEAR=(x1,y1,Z1/x2,y2,z22/x3,y3,z3/etc.....) 


This option gives the value of electric field components in 
the near field for the antenna at the point or points 


specified. 





OUTPUT (FARF=45,50, 25,50) 


9. PLOT This card will produce normalized polar plots in the 
specifed plane for the stated option. The plane is 

___" or "THETA=____". The 
PLOT CARD overrides the limits of the OUTPUT CARD for the 
same option. If only a normalized pattern is of interest, 


specified by stating either "PHI=_ 


only a PLOT CARD is necessary. If a table of values and a 
normalized pattern is desired, both a PLOT CARD and OUTPUT 
CARD must be used. Only one PLOT CARD is permitted per 
data block. The following pattern plots are available: 


FAR FIELD/plane 


This option will plot the far field intensity for each 


component of the electric field. 


BACKSCATTERING/plane 


This option will plot the normalized magnitude of each of 


the elements of the polarization scattering matrix. 
BISTATIC/plane 


This option will plot the normalized magnitude of each of 
the elements of the polarization scattering matrix produced 
by the incident plane wave stated by final limits of the 


backscattering option of the output request. 





PLOT (FARF/THET=90) 





113 


10. GROUND This card 1s used to describe the ground 
parameters if a ground plane is present. If no ground plane 
is present, the structure is assumed to be in free space or 
the homogeneous medium of the EXTERIOR MEDIUM data card. 


Seven options are available and are defined as: 
PERFECT 


This option will analyze the structure over a perfect ground 


plane. 
GOOD 


This option will analyze the structure over a good ground 
plane where the conductivity of the ground is .02 mhos/meter 


and the relative dielectric constant is 30. 
POOR 


This option will analyze the structure over a poor ground 
plane where the conductivity of the ground is .001 


mhos/meter and the relative dielectric constant is 4. 
SEA 


This option will analyze the structure over salt water where 
the conductivty of the water is 4. mhos/meter and the 


relative dielectric constant is 80. 
HEIGHT=value in meters 


This option will analyze the structure with origin of the 
GEOMETRY card this height above the ground plane. The 
lowest point of the structure must not lie below the ground 


plane. It may lie on the ground plane. 
CONDUCTIVITY= value in mhos/meter 


This option is used to state the value of conductivity of 
the ground plane if the default values mentioned above are 


not utilized. 


DIELECTRIC= value 


This option is used to state the relative dielectric 


constant of the ground plane if the default values mentioned 
above are not utilized. 


GROUND (HEIG=10/COND=.002/DIEL=10) 





11. INTERVAL FOR CALCULATION This card is used to state the 


number of intervals to be used for calculating the elements 


of the impedance matrix with Simpson's-rule integration. A 


large value for the number improves the accuracy at the 


expense of greater execution time. For most problems a 


suitable combination of speed and accuracy 
a value of four, the default value. If the 
closed-form impedance expressions in terms 


exponential integrals is desired, set this 


INTERVAL=value 


INTE (6) 


is obtained with 
rigiorous 
of the 


value to zero. 


12. GENERATOR This card is similiar to the FEED CARD in use, 
except that the segment numbers are stated instead of the 
node numbers. This is useful if three or four segments meet 
at a node. The positive terminal of the generator is 
connected to the specified segment such that current is 
forced in the the positive direction. The positive direction 
of current flow is from the first stated node number of that 
segment toward the second stated as ordered on the 
DESCRIPTION CARD. 





GENE(2,-5,-90/4, .5,+90) 
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13. IMPEDANCE This card is similiar to the LOAD CARD in use, 
except that the segment numbers are stated instead of the 
node numbers. AS with the GENERATOR CARD, this is used if 
three or four segments are connected to a node. The 
impedance will be connected to the positive terminal of the 
specified segment. The format of this card is the same as 


the LOAD CARD. 





IMPE(1,120,-45/3, 120, +45) 
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14. CHANGE This card at the end of the data block signals 
the program that the following data cards are changes to the 
previously read data, for the next run. If a "CHANGE CARD" 
is used, the outputs must be requested again in the next 


data block. 


15. END This card signals the program that this is the end 
of a data block series and to reinitialize data for the next 


problem. An “END CARD" cannot be used with a "CHANGE CARD". 


16. STOP This card signals the program that all of the data 
cards have been read and to terminate itself when execution 
is completed. This card must be used as the last card in 
place of the "END CARD" of the last data block series. A 
"STOP CARD" cannot be used with an " END CARD" in the same 
data block. 
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C AN EXAMPLE PROBLEM 

C V ANTENNA 

WIRE (RADIUS 

GEOM(0,-~ a ep ane » 0949 +e 09/09 070/07 02099 009/05 e18y018) 
cate eee yt CG 25) 

FEED USN 

OUTPUT (FARF=45, 50,65, 80/STEP=5) 

CHANGE 

OUT PUT (BIST=455 455 45) 45/BACK=0,0,)10,12) 

OUTPUT (CUR 

CHANGE 

C CHANGE STRUCTURE SHAPE TO DIPOLE 
GEOM(0,—0 25 90/0 9—-012530/09030/05 0125 40/0 502530) 
PLOT (FARF/PHI=90) 

GROUND ( HEIGHT =.25/G000) 


THE ABOVE DATA DE 
FO 


= 
FT OY 


ngs THE OUTPUT ON THE 
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